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Abstract 

We calculate the polarization of radiation from thick accretion disks with vertically averaged 
global magnetic field. The polarization arises as a result of the radiation scattering by free electrons 
in magnetized plasma of a disk. We consider as a basic effect the Faraday rotation of polarization 
plane along the photon propagation in a magnetized disk. The various models of optically thick 
accretion disk with a vertically averaged magnetic field are considered. The main goal of this paper 
is to obtain simple asymptotic formulae for the polarization of radiation when the Faraday rotation 
angle $ > 1 at the Thomson optical length r « 1. The results of our calculation allows us to 
estimate the magnetic field magnitude near the marginally stable orbit region of a black hole via 
the data of polarimetric observations including the expected future X-ray polarimetric observations. 
The wavelength dependence of polarization is strongly dependent on various models of an accretion 
disks and thus allows to choose a real model from the polarization data. 


1 Introduction 

An accretion disk is one of the basic elements for structure of accreting flows around the compact objects 
- neutron stars, black holes, quasars (QSO) and active galactic nuclei (AGN). Due to the absence of the 
axial symmetry with respect to a line of sight, the total integrated radiation emerging from accretion 
disks should be polarized via the electron scattering process. The presence of a magnetic field gives 
a new effect provided the Faraday rotation of a polarization plane along a photon mean free path in 
scattering medium. Namely, a nontrivial wavelength dependence of polarization arises when the Faraday 
rotation angle at the Thomson optical depth r (see Gnedin and Silant’ev, 1997) 

* =04 (A) A)"™" (i) 

is sufficiently large. Here, A is the radiation wavelength and 6 is the angle between the line of sight and 
magnetic field B. 

The key point is a solution of the Milne problem in the case of magnetized atmosphere. This problem 
was considered by many authors, who obtained the numerical solutions (Silant’ev, 1994, 2002, Agol and 
Blaes, 1996, Agol et al., 1998, Shternin et al. 2003). 

The main goal of this paper is to obtain the analytical asymptotic formulae for the polarization of 
radiation and to use them for analysis of various models of an accretion disk with a global magnetic 
field. Here, it is more convenient to use the simple approximate formulae for the Stokes parameters of 
emerging radiation from magnetized, plane-parallel optically thick atmosphere, obtained by Silant’ev, 
2002 . 
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2 The Degree of Polarization from Optically Thick Accretion 
Disks with a Strong Magnetic Field 


Silant’ev, 2002, obtained the quite simple analytical approximate formulae for the degree of polarization 
of radiation emerging from optically thick accretion disk which for considered case vertically averaged 
magnetic field acquires the form 


( x _ 1 - 9 1 ~ M _1_ 

PlW 1 + 3 ' J(/i) ■ ^/(l - kg) 2 + 5V ’ 

where g = cos0, 0 is the angle between the magnetic field directed along the outward normal (z-axis) and 
the radiation emerging from a disk. The function J(/i) describes the angular distribution of emerging 
radiation. This function and the numerical parameters g and k are tabulated by Silant’ev, 2002. The 
parameter of the Faraday depolarization can be introduced from Eq. 0 


6 = 0.8 




(3) 


We consider the case of dominant electron scattering process inside an accretion disk (fc = 0, g = 
0.83255). 

Photons escape the optically thick disk basically from the surface layer with r « 1. If the Faraday 
rotation angle U/ corresponding to this optical length becomes greater than unity, then the emerging 
radiation will be depolarized as a result of the summarizing of radiation fluxes with very different angles 
of Faraday’s rotation. Only for directions that are perpendicular to the vertical magnetic field the 
Faraday rotation angle is too small to yield depolarization effect. Certainly, the diffusion of radiation 
in the inner parts of a disk depolarizes it even in the absence of magnetic field because of multiple 
scattering of photons. The Faraday rotation only increases the depolarization process. It means that 
the polarization of outgoing radiation acquires the peak-like angular dependence with its maximum 
for the direction perpendicular to magnetic field. The sharpness of the peak increases with increasing 
magnetic field magnitude. The main region of allowed angles appears to be ~ 1/5. 

Another the very important feature characterizing the polarized radiation is the wavelength depen¬ 
dence of polarization degree that is strongly different from that for Thomson’s scattering. 

For strong magnetic field magnitude (or large wavelength) when ^ > 1 the simple asymptotic 
formulae takes place 


Pi{B) 


Pi{B = 0) 
5g 


(4) 


where pi(B = 0) is the classical Cliandrasekhar-Sobolev polarization due to pure electron scattering 
without magnetic field. 

Silant’ev, 2002, has considered also the other forms of magnetic field distributions, namely, pure 
radial and toroidal magnetic fields. In both cases the Eqs. 0 and 0 take the same form with 
changing parameter 5g to <5(1 — g 2 ) 1 ^ 2 . For chaotic magnetic field the analytical asymptotic dependence 
does not radically change, i.e. pi(B) ~ 1/5. 

Here we restrict ourselves to the case of an accretion disk with vertically averaged magnetic field. 
Others distributions of a magnetic field will be considered in separate works. 


3 Models of Thick Accretion Disks with Vertically Averaged 
Magnetic Field 

First of all we start with a short review of present models of an accretion disk with vertically averaged 
magnetic field (see, for example, Campbell, 1997, Casse and Keppens, 2002, Li, 2002, Wang et ah, 2003, 
Romanova et ah, 2003, Pariev et ah, 2003, Turner et ah, 2003. 

Campbell, 1997, in his book ’’Magnetohydrodynamics in binary stars” considered the case of a simple 
dipole magnetic field. For radial distances from gravitation centre r z he presented the magnetic field 
distribution as 



p 



zB z 
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The processes of spin-up and spin-down of magnetized stars with accretion disks and outflows are 
considered by Romanova et ah, 2003. They has made three-dimensional simulations of disk accretion 
to an inclined dipole. It was found the radial dependence of vertically averaged magnetic field as 
B z ~ r -5 / 4 . Using their result one can estimate the vertical magnetic field magnitude as 


B z 
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where M is the mass of a compact object (a black hole) and M is an accretion rate magnitude, R g is a 
gravitational radius. 

Li, 2002, obtained for self-graviting magnetically support disks the following radial dependence of 
the vertically averaged magnetic held 




(7) 


where £ = M/2irrv r is the surface density of a disk. If the equipartition condition would be accepted 
it provides the following dependence 


B z 


10 2 


f io 8 M e \ 
\ M J 


(—) 

Y 10 26 g/s J 


1/2 



( 8 ) 


Casse and Keppens, 2002, developed a model depending on the central region of a disk, i.e. on 
/3 = 8nP gas (z = 0 )/B 2 (z = 0). They obtain the following expression for vertically averaged magnetic 


held 
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(9) 


Pariev et al., 2003, extended the Shakura-Sunayev approach to strongly magnetized accretion disk 
model. They assume that the radial dependence of the vertically averaged magnetic held in the disk is 
described by the power law 


B z — B io 



( 10 ) 


where B io is the strength of the magnetic held at 10 R g and 77 > 0 is some constant. They made 
numerical calculations and presented plots of radial structure and emission spectra from the disk in the 
region when it is optically thick for four choices of basic parameters: 77 = 5/4, B 10 = 3 • 10 3 G; 77 = 1, 
B 10 = 3 • 10 3 G; 77 = 1.4, B w = 5 • 10 3 G; 77 = 5/4, B w = 700G and for the mass of supermassive black 
hole M = 10 8 M 0 . 

Liu et al., 2003, developed also the model that appears quite close to Pariev et al. model. Liu et ah, 
2003, obtained the following expression for the averaged magnetic held 


B,. 7 {Jvy V! ° (VkT ^‘ (11) 

Here, a is the Shakura-Sunayev viscosity coefficient and Me is the Eddington accretion rate, Me = 
1.4 • 10 26 (M/10 8 Mq)<7s _1 . One can see that their radial dependence of the averaged magnetic held is 
practically the same as in Pariev et al. model with 77 = 5/4. 

Turner et al., 2003, developed the model of radiation pressure supported accretion disk. In their 
model the vertically averaged magnetic held was presented by the relation 


B z > 10 8 



( 12 ) 


where 77 = 3/4. 

The magnetic held was also included in the standard models of an accretion disk developed by 
Shakura and Sunayev 1973 and by Narayan an Yi, 1995 (ADAF model), without detail specification of 
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its geometry. For these both models the index 77 = 5/4. For example, Shakura and Sunayev gave the 
following expression for the magnetic field strength 


B z — B 10 
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It means that Bio = B\o(M, M ) and 77 = 51/40 ss 5/4. 
The effective temperature of a disk is determined by 
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(14) 


The disk radiates as a black body. It means that the peak spectral wavelength is A = 0.29/T e . The 
last expression allows to get the relation between radial distance from a compact object and the peak 
spectral wavelength corresponding to this distance 
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As a result one obtains 
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The strong depolarization takes place only for those wavelengths which yield 5 > 1. 


(16) 


(17) 


4 The Faraday Depolarization Parameter 5 for Various Disk 
Models 


We calculate now the value of the Faraday depolarization parameter 5 for various models of an accretion 
disk. We are starting with the situation of an accretion disk in close X-ray binaries and consider two 
cases of magnetic field power law: 77 = 3 (dipole field) and more popular case when 77 = 5/4. For 77 = 3 
one obtains from Eq. m 


5 = 0.5- 10 -11 
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(18) 


For optical (V band) photons 5 ~ 2 • 10 11 ( \qm q ) (V/ 10 17 ( 7 /s^ B w . It means that for depo¬ 
larization of optical photons too much magnetic strength B > 10 11 (G') is required. For X-rays with the 
energy E = 1 keV (A = 1.2A) the magnetic field strength B > 10 4 (G) provides the depolarization of 
X-ray photons from optically thick accretion disk. 

For 77 = 5/4 and for V-band of radiation the depolarization parameter becomes equal to 5 = 1.35 • 

/ n5/ 6 / . \ —5/12 

10 -5 ( 1 O m 0 ) (V/10 1 '.g/sJ Bio, and magnetic field strength at the inner radius of a disk B > 

10 5 (G) is quite sufficient for strong depolarization. The result of calculation parameter <5 for other values 
of the index 77 are presented at Table 1 (V-band photons) and Table 2 (X-ray photons with the energy 
IkeV). The most surprising result arises when 77 = 3/2. In this case the depolarization parameter 5 
does not depend on the wavelength (energy) of photons from a disk. It means that the polarization will 
not depend on the wavelength (energy) of photons as it happens in the classical Sobolev-Chandrasekhar 
case. However, the value of polarization will be less due to the depolarization effect. 

One can predict also that, if 77 > 3/2 the hard energy photons will be tested by most strong depo¬ 
larization than soft energy photons, and the polarization spectrum will drop in the hard energy region. 
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On the other way, if 77 < 3/2, the soft energy photons will be stronger depolarized and the spectrum 
will drop in the soft energy region. 


5 The Asymptotic Wavelength Dependence of Polarization 
from the Magnetized Standard Accretion Disk 

Now we present the asymptotic formulae for the polarization degree using the result of analytical cal¬ 
culations by Silant’ev, 2002. (see Eqs. 0 and 0 ). 

According to Eq. 0 the asymptotic wavelength dependence of the optical polarization is determined 

by 
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Thus the asymptotic wavelength behaviour of the polarization allows to derive the value of the index 
of the power-law radial dependence of the vertically averaged magnetic field for the standard magnetized 
accretion disk. The results of the calculations for the various models of the standard magnetized disk 
are presented in the Table 3. 

The basic result from the Table 3 is the radical change of the wavelength dependence of the polar¬ 
ization at the region of the index values 77 = 1.5. For the more strong radial dependence of the magnetic 
field (77 > 3/2) the polarization is increasing with the increase of the wavelength. For the more weak 
radial dependence (77 < 3/2) the polarization is increasing in the short wavelength region, i.e. pi ~ A _ra , 
where n = f(rf). It means that for standard magnetized accretion disk the only fact of the qualitative 
wavelength dependence of polarization (increase or decrease the wavelength) gives an evidence of the 
radial dependence of the magnetic field inside an accretion disk. 


6 Summary and Conclusions 

We presented the formulae for the asymptotic wavelength dependence of the polarization of radiation 
from thick accretion disks with vertically averaged magnetic field. The polarization arises as a result of 
the scattering of radiation by electrons in magnetized plasma of a disk. The effect of the Faraday rotation 
of polarization plane in a magnetized disk is taken into account. We considered the different models of the 
thick accretion disk with the vertically averaged magnetic field. The asymptotic wavelength behaviour 
of the net polarization displays very strongly the power-law radial dependence of the vertical magnetic 
field. It appears that for more strong dependence (the index of the power-law radial dependence 77 > 1.5) 
the hard energy photon will be tested by more strong depolarization than soft photons. For more smooth 
radial dependence of the magnetic field the soft energy photons will be stronger depolarized and spectrum 
of polarized radiation will drop in the soft energy region. This fact displays the qualitative estimation 
the power-law radial dependence of the magnetic field using only the simple observations of the super 
scale dependence of polarization respect to wavelength of photon (increase or decrease depending on the 
wavelength). 
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Table 1. 

The values of the magnetic field strengths B\o for various magnitudes of the power law index of the 
magnetic field distributions in the standard accretion disk for V-band photons, if S > 1. 
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Table 2. 

The values of the magnetic field strengths B io for various magnitudes of the power law index of the 
magnetic field distributions in the standard accretion disk for X-ray photons ( E x ss 1 keV), if S > 1. 
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Table 3. 

The asymptotic waveleghth dependence of the polarization for the standard magnetized accretion disk. 
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